Hysteresis in an Ising model with mobile bonds 
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Hysteresis is studied in a disordered Ising model in which diffusion of antiferromagnetic bonds 
is allowed in addition to spin flips. Saturation behavior changes to a figure-eight loop when dif- 
fusion is introduced. The upper and lower flelds delimiting the figure-eight are determined by the 
Hamiltonian, while its surface and the crossing point depend on the temperature and details of 
the dynamics. The main avalanche is associated with the disappearance of hidden order. Some 
experimental observations of figure-eight anomalies are discussed. It is argued they are a signal of a 
transient rearrangement of domain couplings, characteristic of amorphous and/or magnetically soft 
samples, and similar to evolution of kinetic glasses. 

I. INTRODUCTION 

Magnetic hysteresis is a classic example of metastability at experimentally easily accessible time scales. Its ready 
appearance under simple conditions makes it not only technologically important, but also a theoretical challenge. 
In addition to macroscopic approaches, relying on multiply-valued equations of state, there has been considerable 
work done over the past decade in Monte Carlo simulations of lattice systems. The high reproducibility of hysteresis 
curves, and the fact that among their many observed shapes most are fairly simple, imply that hysteretic delay involves 
processes similar to ordinary thermalization, which can be modeled by a simple random walk. It is thus reasonable to 
ask, which features need to be added to the walker in order to reproduce some particular aspect of the experimental 
observations. 

The Ising model, in the widest sense, has evolved as an important test-bed for this kind of question. Its spins are 
naturally associated with magnetic domains, and can be coupled by different physically motivated interactions. It 
turns out that random fields with cluster updates give rise to Barkhausen noise |lj], local spin inversion symmetry is 
needed for memory effects 0, random couplings produce multicycles at low temperatures 13 , while long-range and 
anisotropic interactions give domain patterns very similar to those in some real thin films U • 

The present work reports one more result of this kind: allow local diffusion of antiferromagnetic (AF) couplings, 
and the saturation part of the hysteresis curve will change, from the usual horizontal line to a figure-eight form. Such 
anomalies are sometimes observed in experiments on thin films and tapes, when the coercive fields are of the order 
of 100 Oe. We are aware of only one paper, however, where the authors have specifically commented upon them j5|. 
It is known from a diffraction experiment, albeit with domain structure carefully controlled by micromachining, that 
spatial evolution of AF couplings can occur under driving by a field We suggest here that a physically similar, 
but random, mechanism may be responsible for figure-eight shapes in saturation curves, in the usual situation, where 
sample domain structure is left to chance. 

II. THE MODEL 

We use the same model as before 0- Take the two-dimensional short-range bond-disordered Ising model, or 
±J-Edwards- Anderson (EA) model 

H = -Y. J^,S,S,-BY,S^ (1) 

<i,j> i 

where Jij = ±J and Si ~ ±1. The model is subject to Glauber |2J dynamics: spin is flipped at random, and the move 
is accepted if the criterion 



1 + exp(A/T) 



> c 



is satisfied, where A is the energy balance of the trial move, T the temperature, and c a uniformly distributed random 
number, < c < 1. Between each two spin trials there is a bond trial: a positive and negative bond impinging on 
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FIG. 1: Hysteresis in the model Q, with (full curve) and without (dotted curve) bond updates. Time spent at each field is 
ten updates per site. The initial parts of the loops have been erased for clarity. Note that the half-width of the main loop is 
the same as the range of each figure eight, namely 2 J. Open circles: see Fig.|3] 



the same site, and chosen at random, are allowed to exchange places, subject to the same criterion as the spin trials. 
The concentration paf of antiferromagnetic (AF) bonds is taken to be 50%, unless otherwise noted. 

This dynamics causes a significant annealing of the sample, and after a transient of 100 updates per site it 
enters a long-lived state characterized by glass-like dynamics, with a second, much longer relaxation time [J. The 
'glassiness' of the metastable state is not topological: the underlying bond distribution, viewed at any instant in time, 
has low frustration, so much so that it can be mapped onto a disordered ferromagnet with a finite transition 
temperature However, when the bonds move, the spins are prevented from settling down into any given ordered 
arrangement. In this sense the situation is more reminiscent of the 'geometric' frustration of vitreous liquids, than 
of spin glasses: configurations are locally relaxed, but the spin correlation decays over time even at low temperature. 
As long as the number of F and AF bonds is kept equal, the equilibrium zero-field magnetization is zero. 

The equivalent disordered ferromagnet is obtained by local gauge transformations which change the sign of bonds 
whenever the plaquette is not frustrated. The disorder in the ferromagnet comes from the residual AF bonds, related 
to the frustrated plaquettes, which are relatively few, so a net magnetization can appear, revealing the formerly 
hidden order. (This is different from a dilute ferromagnet, where some bonds have been assigned J = 0.) The degree 
of frustration, hence the order associated with the hidden ferromagnet, can change as a function of external field, 
because of bond diffusion. This is specific to the present model. In the usual situation without bond updates, the 
frustration is fixed for a given sample, and when paf = 50%, the equivalent ferromagnet has no phase transition at 
finite temperature, and neither has the spin glass. 



III. MAIN RESULT 



The model is subjected to a periodic external field, starting with _B = and increasing to B = 4J, then reducing 
the field to B = — 4J, and finally closing the curve. At each point in the field, we take 10 updates per site (meaning 
10 spin and 10 bond updates in alternation). The dependence on this number (sweep rate) will be discussed in 
the next section. Two field steps were used, AS ~ 0.04J and AS ~ 0.08J, with no discernible difference in the 
results. We note that for our two-dimensional square lattice, i? = 4J is the limit at which saturation is enforced, 
at a magnetization determined only by the temperature. The outcome is shown in Fig. ^ compared to the case 
without bond updates. The latter has been studied before fo'^ the 3D case. Its hysteresis curve consists of a 
set of box-like regions, with constant magnetization alternating with large avalanches. This is easily understood, 
because the system is essentially in the clean limit. Namely, an individual Ising spin has only a finite, discrete set 
of energy values available. As soon as the external field passes through one of these values, all the respective spins 
are destabilized at the same time, and there is nothing to stop the ensuing avalanche. This result corresponds to low 
temperature (T = 0.1 J), on which we concentrate in the greater part of the article. Note that with a discrete set 
of bond strengths, the system is gapped, so the temperature T = O.IJ gives nearly the same results as T = O.OIJ, 
except that the figure eight is slightly shifted, an insignificant effect to the purposes of this article. 

Bond updates are now introduced, alternating equally with spin updates. It is no longer necessary to average 



3 




FIG. 2: Magnified figure-eight anomaly of Fig.0 fitted (tliick curves) to with a — —0.07. The value of a is the same in the 
ascending and descending parts. Circles are raw data over four cycles. Inset: dependence of the exponent on the temperature. 

over samples, because bond movement makes the whole sample space available to a single walker, or in technical 
terms, it induces self-averaging. This is so efficient that Fig. ^ simply gives the raw data for five cycles, showing how 
the corresponding points fall on the same curve, meaning the error bars are smaller than the line thickness. The 
simulations were mainly on 100 x 100 lattices, with spot-checks to make sure that 50 x 50 and 200 x 200 gave the 
same results. The curves were also independent of the random initial conditions, given that field cycling began after 
the above-mentioned fast transient (50-100 MCS per site) died out. 

Bond movement introduces two new features in Fig.^ First, the two inner boxes (around the origin) open up into 
a single large one, with avalanches at -B = ±2 J. This is similar to what happens when a random static site field is 
added to the model without bond movement, or for a Gaussian distribution of bond strengths. It can be explained 
by the disappearance of spins intrinsically in an indifferent state, surrounded by the same number of positive and 
negative bonds, which are responsible for the avalanches at zero field. However, the reason for the disappearance 
is somewhat different in the two cases. Introducing site disorder of some characteristic magnitude K will make the 
formerly indifferent spins unstable at some new field values B « ±K ^ 0, leaving few to participate in the avalanche 
at _B = 0. In the present model, on the other hand, bond movement induces significant relaxation of the system, 
meaning most spins are in a locally minimal configuration, so again the number of indifferent ones is low. 

The second new feature is somewhat unexpected. The former small box at 2J < B < 4J is now a distinct figure- 
eight loop. Both the forward (increasing B) and backward (decreasing B) parts of the magnetization curve have the 
form of relaxation curves. They are given by power laws in B with small exponents, as shown in Fig. [21 By contrast, 
the temporal autocorrelation of the magnetization is a stretched exponential, as in other similar models. Hence the 
M{B) curves cannot be interpreted as a simple temporal evolution, trivially due to the fact that B changes before the 
system has had time to equilibrate. The bond updates introduce a new relaxation mechanism, with a scale comparable 
to the scale of the main cycle. The exponent of B becomes (absolutely) smaller as the temperature increases, which 
reflects the figure-eight becoming thinner and flatter. The temperature dependence of the exponent is given in the 
inset to the figure. 

The main avalanche in this model is microscopically different than in models with quenched disorder, because in 
the presence of bond updates, the frustration of the lattice can vary as a function of the field. This is shown in 
Figure 13 In the region of the main loop, the equivalent disordered ferromagnet has a finite transition temperature. 
The concentration of frustrated plaquettes rises sharply with the main avalanche, and passes through the critical 
value, above which the ferromagnet has Tc — 0. In this way the spin avalanche is connected with a ground-state 
transition in the bond lattice. For B < 2 J, the instantaneous bond configuration can support an ordered spin state 
at finite temperature; above this value, it cannot. 

In Fig. 21 sample configurations are shown immediately before and after the main avalanche. In the ordered state, 
there are islands of spins with the same orientation (dot: up, circle: down), connected by ferromagnetic bonds. In the 
disordered state, the islands of down-spins (circles) have broken up, and each is preferentially connected with up-spins 
by 'crosses' of AF bonds. 
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FIG. 3: Left: concentration of frustrated plaquettes as a function of the field, for the full curve in Fig. Horizontal line: 
critical value pc = 0.31 at which Tc = for the equivalent disordered ferromagnet. Right: transition temperature of the latter, 
obtained from the values in the left panel according to Ref. flU . 
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FIG. 4: Configurations before and after the main avalanche. Left: ordered state, lower open circle in Fig.Q Only ferromagnetic 
bonds are shown. Right: disordered state, upper open circle in Fig. Only AF bonds are shown. Circles: down-spins. Dots: 
up-spins. 



IV. PHENOMENOLOGY OF THE MODEL 



In this section, we show how the basic curve in Fig. ^ changes with various aspects of the model dynamics. There 
are three to which it should be at least moderately robust: temperature, ratio of spin to bond updates, and sweep 
rate. 

In Figure 13 we observe that the figure-eight evolves with temperature on the scale of J. At T = 0.5 J, the first 
maximum of the figure eight has dropped below the reverse curve, so there is now an open narrow loop at i? > 2 J, 
which is still wider at the middle than at either end. The narrowing of the main part —2J<B<2Jis significantly 
milder than in the case without bond updates jl^, so that, overall, their introduction makes the hysteresis less 
sensitive to temperature increase. 

In Figureini we show that the figure-eight feature is not a transient effect due to too fast cycling. It visibly persists 
to up to 500 updates/site at each field, although it becomes thinner. We know that in the model's two-time dynamics, 
the transition between the fast (transient) and slow (metastable equilibrium) component occurs at 50-100 updates 
per site, so the feature survives into the intrinsic 'glassy' regime. Observed figure-eight anomalies are indeed thinner 
(relative to the main loop) than suggested by Fig. ^ This agrees with the idea that experimentally, magnetizations 
at each point of the loop are measured after any fast transients have died out. We find, nevertheless, that the 
power-law relaxation of M with B is found in all three curves in Fig. ^ so it is not wrong to use fast cycling for 
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FIG. 5: Temperature dependence of the hysteresis. Thick full line: T — O.IJ (Fig. 0. Dotted line: T = 0.3J. Broken line: 
T = 0.5J. Thin full line: T = J. 




FIG. 6: Dependence of the figure-eight feature on the time spent at each field. Full line: 10 updates/site (Fig.Ql. Dotted line: 
100 updates. Broken line: 500 updates. The field step is AB = 0.04J in all three curves. 



model investigations. Only the numerical value of the relaxation exponent changes, decreasing with slower cychng. 
Increasing the number of updates at each field also gets rid of the steps at i? = in Fig. without changing the 
shape and size of the main loop. 

In Figured we show perhaps the most unexpected outcome of the phcnomenological investigations. We vary the 
relative frequency of spin and bond updates, guided by the idea that they are not equally common in real materials. 
Based on past modelling experience, one would assume that bond updates should be rarer, so it matters that the 
figure-eight is not too sensitive to their rate. Indeed it is not, but it also turns out that the system responds much 
the same whether it is the rate of spin or bond updates that is decreased. Hysteresis curves in which they differ 
by a factor of ten either way fall practically on top of each other. Bond updates are therefore efficient in mutually 
equilibrating the spin and bond subsystems. They are out of equilibrium with respect to the asymptotic state, but 
not with respect to each other. Figure [7| also shows a dependence on the total number of updates at each field value, 
i.e. the sweep rate, like in Fig.lHI The curves with opposite ratios of spin to bond updates are similar only if the sum 
of the updates is the same. The thin full line has the same ratio of the two as the dash-dotted line, but shows the 
effect of faster cycling, with a more pronounced figure eight, and noisier magnetization. (Note that the iV = 10 curve 
in Fig. El has the same number of updates as the thick full fine in Fig.d) 

In fact, it turns out that the ±J-Edwards- Anderson model is extremely sensitive to the introduction of bond 
updates. Take the T/J — 0.3 curve in Fig.jS] obtained with a 1:1 ratio of bond to spin updates. Going down to 1:200 
barely changes it: it becomes something like the T/J = 0.5 curve in the same figure. But turning off bond updates 
entirely collapses the main loop to the two boxes shown in Fig. ^ 

Finally, we find that the figure-eight is quite sensitive to noise in the bond couplings. A Gaussian spread with 
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FIG. 7: Dependence of the figure-eight feature on the ratio of spin to bond updates, at T = 0.05J. Raw data for five cycles 
are shown for each curve. See the text for details. 



CT — 0.05 J centered at ± J is sufficient to flatten it almost completely. Thus while it can survive some noise in the 
bond distribution — i.e. is not a singular effect — the figure-eight anomaly should at present be regarded as specific 
to the ± J model. 



V. EXPERIMENTAL EVIDENCE 



Experimental observations of the figure-eight anomaly are not very rare, but it is rare that anyone comments on 
them; in fact we found only one paper in which this was explicitly mentioned as an unexplained phenomenon j^] . This 
may partly be because such forms can be dismissed as artefacts, due to the integration of a noisy differential signal 
dM / dH in the saturation region, where the derivative is small. Another possible reason is that they are sometimes 
observed 13] along with a number of other 'misfeatures,' such as a marked absence of inversion symmetry relative to 
the origin in the {B, M) plane, so it is not clear that they should be singled out for attention. 

In Ref. jlJl, annealed amorphous samples of Co66Fe4Bi5Sii5 with a recrystallized surface were studied by the 
magneto-optical effect. The surface was chemically etched to expose the amorphous interior of annealed samples, and 
the changes in the hysteresis curves recorded. After etching for ten minutes, the magnetic response became softer and 
a thin but distinct figure-eight appeared in the saturation region. Since the laser is a local probe, it was also possible 
to obtain hysteresis curves at different physical positions on the as-quenched sample surface. Some of these also show 
a figure eight, which the authors do not comment. The width of the measured 'boxy' main loop (large avalanche!) 
in the field was of the same order as the range of the figure eight, just as found in the model calculations here (see 
Fig.P). 

In Ref. 0, a thin-film grating was subject to a perodic field and the magnetic order observed by X-ray scattering. 
A perfectly reproducible AF ordering was observed at both values of the coercive field, i.e. as the sample passed 
through zero magnetization. This does not mean that the positions of the effective AF (dipolar) couplings among the 
domains has been changing, in fact quite the opposite, it is proof that it is fixed. The interesting observation from our 
point of view occurred during sample preparation. Direct observation of magnetic domains in 'fresh' samples showed 
they did not follow the lines of the grating, but were in fact quite random, so that the grating geometry could barely 
be discerned. A careful demagnetization process, with a slowly vanishing field, was necessary to arrange the domains 
along the grating lines. During this process, the effective couplings were evidently rearranged in space. Once they 'fell 
into line,' they did not evolve any further. We suggest a similar evolution is happening in ordinary 'fresh' samples, 
without control of domain positions. First there is some rearrangement of domains under field cycling, amounting to 
a spatial evolution (diffusion) of AF bonds. As the sample ages, the domain positions should become fixed, pinned 
perhaps by impurities, and the figure-eight anomaly should disappear. This scenario is corroborated in Ref. |l3j . 
whose authors observed a hardening of the magnetic response and disappearance of the anomaly after eight hours of 
annealing, which recrystallized the surface. When the sample was etched as described above, the anomaly reappeared. 
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VI. DISCUSSION 

Hysteretic effects are to be expected whenever the approach to thermal equiUbrium is somehow hindered. One 
widely studied way in which this can happen is when the equilibrium state is itself disordered, and the system 
takes a long time to find this state among the many suboptimal, similarly disordered ones. This is the case of spin 
glasses, especially in two dimensions where the spin-glass transition temperature is zero. A different way to hinder 
equilibration is provided by kinetic glasses, where the underlying equilibrium state is ordered, but is never reached 
because of kinetic obstruction. 

The use of the disordered Ising model with fixed bond distribution to study hysteresis assumes that the metastability 
involved in magnetic hysteresis may have something in common with the metastability of spin glasses. This is tacitly 
assumed in all such investigations of hysteresis with fixed bond distributions. The present work investigates the 
alternative possibility, that the delayed thermalization in magnetic hysteresis may be more like the metastability 
of kinetic glasses. The figure-eight anomaly was identified as a distinct and sensitive signal of this latter scenario. 
Experimental observations Q allow that a real-space redistribution of effective AF interactions between magnetic 
domains can be driven by an external field, with response times within the experimental window. One would intuitively 
expect this rearrangement capability to decrease as the sample ages. Indeed this seems to be the case in the experiment 
mentioned above 

The above observations open some interesting theoretical issues, as well. One is a possible connection between the 
hidden order-disorder transition in the bond background at the main avalanche, and the appearance of the figure-eight 
anomaly in the state disordered by the field. A related issue is to understand the origin of the power-law relaxation 
of the magnetization as a function of the field. Another is the dependence of the anomaly on the distribution of bond 
couplings — we know, for example, that there is no anomaly for a unimodular Gaussian distribution. These questions 
are currently under investigation. 

To conclude, we have given a simple interpretation of the figure-eight anomaly, sometimes observed in magnetic 
thin films and tapes in the saturation region, when the coercive fields involved are relatively weak, of the order of 
100 Oe. We believe it is related to a specific transient process, which may be represented as a spatial rearrangement 
of AF (dipolar) couplings in the sample, driven by the field. 
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